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Abstract 


There  are  limitations  and  shortcomings  to  the  Fourier  transform  method  to  detect  accelerat¬ 
ing  targets  because  of  the  phenomenon  known  as  Doppler  smearing.  In  using  a  Pulse  Doppler 
Radar  to  detect  a  non-accelerating  target  in  additive  white  Gaussian  noise  and  to  estimate 
its  radial  velocity,  the  Fourier  method  provides  an  output  signal-to-noise  ratio  (SNR)  that 
increases  linearly  with  the  number  of  pulses.  When  the  target  is  accelerating,  the  Fourier 
method  may  still  be  used  to  detect  the  target  and  estimate  its  median  velocity,  provided  the 
acceleration  is  small  enough  in  the  sense  described  in  the  paper.  For  a  given  acceleration, 
when  the  number  of  pulses  is  increased,  the  output  SNR  of  the  Fourier  method  varies  as 
a  concave  function,  increasing  to  a  maximum  and  then  decreasing,  before  the  method  fails 
catastrophically.  Thus  the  number  of  pulses  and  the  acceleration  have  to  be  matched  to 
achieve  optimum  performance.  Empirical  formulae  for  the  dependence  of  the  optimum  SNR 
and  the  optimum  number  of  pulses  on  the  acceleration  are  given.  The  results  are  shown  to 
be  relevant  to  the  design  of  Generalized  Likelihood  Ratio  Test  (GLRT)  based  detectors  that 
apply  a  search  over  a  grid. 


Resume 


Us  y  des  Smites  et  des  problemes  associes  avec  l’utilisation  de  la  methode  de  transformdee 
de  Fourier  pour  la  detection  des  cibles  qui  accelerent  a  cause  du  phenomene  d’elargissement 
Doppler.  Dans  Putilisation  des  radars  Doppler  a  impulsions  pour  la  detection  de  cibles  non 
accelerante  en  presence  de  bruit  blanc  et  l’estimation  de  la  velocite  radiale,  la  methode  de 
transformee  de  Fourier  fournit  un  rapport  de  signal  au  bruit  (SNR)  qui  augmente  lineairement 
avec  le  nombre  de  pulses.  Quand  la  cible  est  en  train  d’accelerer  la  methode  de  Fourier  peut 
etre  encore  utilisee  pour  detecter  la  cible  et  sa  vitesse  mediane  en  autant  que  Pacceleration 
soit  suffisamment  petite  comme  montree  dans  ce  rapport.  Pour  une  acceleration  donnee, 
quand  le  nombre  de  pulse  augmente,  le  SNR  resultant  de  la  transformee  de  Fourier  varie 
comme  une  fonction  concave,  augmentant  a  un  maximum  et  apres  decroissant,  avant  que  la 
methode  echoue  de  fagon  catastrophique.  Ainsi  le  nombre  de  pulse  et  Pacceleration  de  la 
cible  doivent  etre  relier  ensemble  pour  accomplir  une  performance  optimum.  Des  formules 
empiriques  donnant  le  SNR  optimum  en  fonction  du  nombre  de  pulse  et  Pacceleration  sont 
donnees.  Les  resultats  sont  relevant  au  developpement  de  detecteurs  utilisant  un  test  de 
rapport  de  probabilite  generalise  qui  applique  une  recherche  a  travers  une  grille. 
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Executive  Summary 


The  Fourier  Transform  is  at  the  heart  of  a  wide  range  of  techniques  that  are  used  in  HF 
radar  data  analysis  and  processing.  Mapping  the  data  into  the  temporal  frequency  domain  is 
an  effective  way  of  recording  the  data  such  that  their  global  characteristics  can  be  assessed. 
However,  the  change  of  frequency  content  with  time  is  one  of  the  main  features  we  observe  in 
HF  radar  data.  Because  of  this  change  of  frequency  content  with  time,  radar  signals  belong 
to  the  class  of  non-stationary  signals. 

One  of  the  central  problems  in  High  Frequency  (HF)  radar  data  is  the  analysis  of  a  time  series. 
The  Fourier  transform  method,  or  Doppler  processing  method,  has  been  generally  used  in 
HF  radar  to  detect  targets  that  are  moving  with  constant  radial  acceleration.  Examples 
of  accelerating  targets  are  manoeuvring  aircrafts  and  missiles.  In  this  report  we  show  that 
there  are  limitations  and  shortcomings  to  the  Fourier  transform  method  to  detect  accelerating 
targets  because  of  the  phenomenon  known  as  Doppler  smearing.  We  show  that  when  the 
target  is  constantly  accelerating,  the  Fourier  method  may  still  be  used  to  detect  target  and 
estimate  its  median  velocity,  provided  the  acceleration  is  small  enough  in  the  sense  to  be 
described  in  this  report.  It  is  shown  that  for  a  given  acceleration,  the  number  of  pulses  cannot 
be  increased  indefinitely  without  resulting  in  catastrophic  failure  of  the  method.  Conversely, 
for  a  given  number  of  pulses,  the  acceleration  cannot  be  arbitrarily  large  without  resulting 
in  catastrophic  failure  of  the  method.  Thus  the  number  of  pulses  and  the  acceleration  have 
to  be  matched  to  achieve  optimum  performance. 


Author.  Thayananthan,  Thayaparan,  Limitations  and  strengths  of  the  Fourier  transform  method  to  detect 
accelerating  targets,  Defence  Research  Establishment  Ottawa,  DREO  TM  2000-078,  November  2000. 
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Sommaire 


La  transforms  de  Fourier  est  au  coeur  d’un  large  eventail  de  techniques  qui  sont,  utilisees 
dans  1  analyse  et  le  traitement  des  donnees  des  radars  a  hautes  frequences.  Le  mappage 
des  donnees  dans  le  domaine  temporel  de  frequence  est  une  fa§on  effective  d’enregistrer  les 
donnees  tel  que  leurs  caracteristiques  globales  peuvent  etre  evaluer.  Cependant,  le  change- 
ment  des  composantes  spectrales  avec  le  temps  est  une  des  caracteristiques  que  nous  obser- 
vons  dans  les  donnees  radar  a  haute  frequence.  Ces  signaux  radars  appartiennent  a  la  classe 

des  signaux  non  stationnaires  a  cause  de  ce  changement  de  composantes  spectrales  avec  le 
temps. 

Un  des  problemes  centraux  dans  les  donnees  des  radars  a  haute  frequence  est  l’analyses 
des  series  temporelles.  La  methode  des  transforme  de  Fourier  ou  la  methode  de  traitement 
Doppler  a  ete  generalement  utilisee  dans  les  radars  a  haute  frequence  pour  detecter  des  cibles 
qui  se  deplace  avec  une  acceleration  radiale  constante.  Exemples  de  cible  qui  accelerent  sont 
des  avions  et  des  missiles  qui  manoeuvrent.  Dans  ce  rapport,  nous  montrons  qu’ils  y  des 
limites  et  des  problemes  associes  avec  l’utilisation  de  la  methode  de  transformee  de  Fourier 
pour  la  detection  des  cibles  qui  accelerent  a  cause  du  phenomene  d’elargissement  Doppler. 
Nous  montrons  que  quand  la  cible  est  en  train  d’accelerer  la  methode  de  Fourier  peut,  etre 
encore  utilisee  pour  detecter  la  cible  et  sa  vitesse  mediane  en  autant  que  l’accelerat.ion  soit. 
suffisamment  petite  comme  montree  dans  ce  rapport.  Ce  rapport  montre  que  pour  une 
acceleration  donne,  le  nombre  de  pulse  ne  peut  etre  augmenter  indefiniment  sans  que  la 
methode  echoue  de  faQon  catastrophique.  D  un  autre  cotee,  pour  un  nombre  donne  de 
pulses,  1’ acceleration  ne  peut  etre  infiniment.  grand  sans  que  la  methode  echoue  de  fa<;:on 
catastrophique.  Ainsi  le  nombre  de  pulse  et  l’acceleration  doivent  etre  relier  ensemble  pour 
accomplir  une  performance  optimum. 


Auther,  Thayananthan,  Thayaparan,  Limitations  et  points  forts  de  la  methode  de  transformee  de  Fourier 
pour  la  detection  de  cibles  qui  accelere,  Center  de  Recherches  pour  la  Defense,  Ottawa,  DREO  TM 
2000-078,  Novembre  2000. 
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1  Introduction 


The  Fourier  Transform  is  at  the  heart  of  a  wide  range  of  techniques  that  are  used  in  HF 
radar  data  analysis  and  processing.  Mapping  the  data  into  the  temporal  frequency  domain  is 
an  effective  way  of  recording  the  data  such  that  their  global  characteristics  can  be  assessed. 
However,  the  change  of  frequency  content  with  time  is  one  of  the  main  features  we  observe 
in  HF  radar  data.  Because  of  this  change  of  frequency  content  with  time,  radar  signals 
belong  to  the  class  of  non-stationary  signals.  One  of  the  central  problems  in  High  Frequency 
(HF)  radar  data  is  the  analysis  of  a  time  series.  The  Fourier  transform  method,  or  Doppler 
processing  method,  has  been  generally  used  in  HF  radar  to  detect  targets  that  are  moving 
with  constant  radial  acceleration.  Examples  of  accelerating  targets  are  manoeuvring  aircrafts 
and  missiles. 

In  this  paper,  the  radar  is  assumed  to  be  of  the  Pulse-Doppler  type,  that  is,  it  sends  out 
a  uniform  train  of  radio  frequency  pulses  and  phase-coherently  receives  their  returns.  It  is 
also  assumed  that  the  receiver  has  pulse  compression  capability  so  that  a  pulse  return  from 
an  isolated  target  can  be  represented  by  a  single  sample  of  the  compressed  pulse. 

A  moving  target  with  low  per-pulse  signal-to-noise  ratio  (SNR)  may  be  efficiently  detected 
by  processing  the  sequence  of  Pulse-Doppler  radar  returns  by  coherent  integration  [  Whalen, 
1971;  Meyer  and  Mayer,  1973;  Wehner,  1987;  Scheer,  1993].  Under  the  simplifying  assump¬ 
tions  of  zero  range-walk  and  constant  radial  velocity,  coherent  integration  is  achieved  by 
taking  the  Fourier  power  spectrum  of  the  sequence  of  complex  amplitudes  of  the  received 
pulses  and  locating  its  maximum.  If  the  maximum  value  exceeds  a  threshold  then  the  pres¬ 
ence  of  a  target  is  declared,  and  the  location  of  the  maximum  is  taken  to  be  an  estimate 
of  the  radial  velocity.  In  an  additive  white  noise  background,  the  Fourier  method  gives  an 
output  SNR  that  is  proportional  to  the  number  of  pulses  [Meyer  and  Mayer,  1973]. 

In  this  paper,  we  relax  the  assumption  of  constant  radial  velocity  and  consider  using  the 
Fourier  method  to  detect  targets  that  are  moving  with  constant  radial  acceleration.  Thus 
the  aim  of  the  paper  is  to  study  the  phenomenon  usually  referred  to  as  Doppler  smearing 
and  quantify  the  performance  of  the  Fourier  method.  In  this  paper,  the  term  acceleration 
may  be  interpreted  as  the  magnitude  of  the  acceleration,  as  the  conclusions  depend  only  on 
the  magnitude  of  the  acceleration.  It  is  empirically  shown  that  when  the  target  is  constantly 
accelerating,  the  Fourier  method  may  still  be  used  to  detect  the  target  and  estimate  its 
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median  velocity,  provided  the  acceleration  is  small  enough  in  the  sense  to  be  described  in 
this  paper. 


It  is  shown  that  for  a  given  acceleration,  the  number  of  pulses  cannot  be  increased  indefinitely 
without  resulting  in  catastrophic  failure  of  the  method.  Conversely,  for  a  given  number  of 
pulses,  the  acceleration  cannot  be  arbitrarily  large  without  resulting  in  catastrophic  failure 
of  the  method.  Moreover,  for  a  given  acceleration,  when  the  number  of  pulses  is  increased, 
the  output  SNR  varies  as  a  concave  function  of  the  number  of  pulses,  achieving  a  maximum 
well  before  the  failure  occurs.  For  a  given  number  of  pulses,  the  ouput  SNR  monotonically 
decreases  with  the  acceleration,  until  the  failure  occcurs.  Thus  the  number  of  pulses  and  the 
acceleration  have  to  be  matched  to  achieve  optimum  performance. 

The  results  are  presented  in  terms  of  a  normalized  value  of  the  acceleration  that  takes 
into  account  the  radar  carrier  frequency  and  the  pulse  repetition  freqency.  This  makes  the 
results  widely  useful.  The  results  may  be  used  to  determine  whether,  in  a  given  scenario 
of  an  accelerating  target,  the  Fourier  method  can  be  used  or  a  more  sophisticated  method 
such  as  the  Generalized  Likelihood  Ratio  Test  (GLRT)  is  required.  In  scenarios  where  the 
Fourier  method  is  not  desirable,  and  where  the  GLRT  method  is  considered,  it  is  shown  how 
the  results  may  be  used  to  determine  the  critical  spacing  of  the  grid  over  which  the  search 
is  applied. 
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2  Signal  Model 


Suppose  the  Pulse-Doppler  radar  sends  out  N  pulses,  one  every  T  seconds,  and  there  is  a 
target  moving  with  constant  radial  acceleration.  Assume  that  the  change  in  range  during  the 
observation  period  of  NT  seconds,  known  as  range  walk,  is  negligible  compared  to  the  radar 
range  resolution  as  determined  by  the  width  of  the  compressed  pulse.  Then  the  complex 
amplitudes  of  the  N  range-compressed  pulses  taken  at  the  range  of  the  target  have  the  form 


r(n)  =  s(n)  +  v(n) 


(1) 


where 


s(n)  =  aej(bo+bin+^b2n2)  (2) 

for  n  =  0, 1, 2,  3, . . . ,  ( N  —  1),  is  the  noise-free  signal,  and  v(n)  is  a  sequence  of  independent 
samples  of  complex  Gaussian  noise  with  mean  zero  and  variance  a2. 

The  signal  parameters  a  and  &o  are  the  target  signal  amplitude  and  phase  respectively.  The 
signal  parameters  b\  and  are  the  normalized  initial  radial  velocity  and  the  normalized 
radial  acceleration  respectively. 


2.1  Normalized  Initial  Radial  Velocity  and  Radial  Acceleration 

The  normalized  initial  radial  velocity  b\  is  defined  as 

“-(f) 

where  u  is  the  initial  radial  velocity  in  meters/sec  towards  the  radar,  T  is  pulse  repetition 
interval  in  secs,  and  A  is  the  carrier  wavelength  in  meters.  Similarly,  the  normalized  radial 
acceleration  b<2  is  defined  as 


(3) 


&2  =  / 


(4) 
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where  /  is  the  radial  acceleration  in  meters/sec/sec  towards  the  radar.  It  can  be  seen  that 
both  of  the  above  normalized  quantities  are  non-dimensional. 

To  obtain  these  relations  note  that  the  decrement  in  range  d,  measured  from  the  beginning  of 
the  observation  interval,  as  a  function  of  continuous-time  t,  is  d(t)  —  ut  +  \ft\  and  therefore 
the  phase  increment  of  a  pulse  as  a  function  of  discrete-time,  or  pulse  index,  is  given  by 


27T 


(5) 


In  the  rest  of  the  discussion,  the  terms  velocity  and  acceleration  refer  to  normalized  radial 
velocity  and  normalized  radial  acceleration  respectively. 


2.2  Per-Pulse  Signal-to-Noise  Ratio 

The  Per-Pulse  Signal-to-Noise  Ratio  is  defined  as 

SNIW  =  (;)2. 


(6) 
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3  The  performance  of  the  Fourier  method  in  detecting 
a  non-accelerating  target 


Although  a  given  signal  can  be  represented  in  many  different  ways,  the  most  important  are 
the  time  and  frequency  representations.  The  majority  of  signals  encountered  in  our  everyday 
life  are  directly  related  to  time.  The  frequency  representations,  on  the  other  hand,  were  not 
popular  until  the  early  19th  century  when  Fourier  first  proposed  the  harmonic  trigonometric 
series.  Since  then,  the  frequency  representation  has  become  one  of  the  most  powerful  and 
standard  tools  for  studying  signals.  By  using  frequency  representations,  we  could  better 
understand  many  physical  phenomenon  and  accomplish  many  things  that  cannot  achieved 
based  on  time  representations.  While  the  time  domain  function  indicates  how  a  signal’s 
amplitude  changes  with  time,  the  frequency  domain  function  tells  how  often  such  changes 
take  place.  The  bridge  between  time  and  frequency  is  the  Fourier  Transform. 

The  signal  is  expanded  in  terms  of  sinusoids  of  different  frequencies  [ Brigham ,  1974;  Bloom¬ 
field,  1976;  Papoulis  1977;  Bracewell,  1978] 

s(t)=-j=J  (7) 

The  waveform  is  made  up  of  the  addition  (linear  superposition)  of  the  simple  waveforms, 
eju}t,  each  characterized  by  the  frequency,  u>,  and  contributing  a  relative  amount  indicated 
by  the  coefficient,  S( of).  S(uj)  is  obtained  from  the  signal  by 

S(w)  =  -^=  j  sit)  e~ju,t  dt  (8) 

and  is  called  the  spectrum  or  the  Fourier  Transform.  The  square  of  the  Fourier  Transform 
|  S(u)  |2  is  called  the  Fourier  power  spectrum ,  which  indicates  how  the  signal  energy  is 
distributed  in  the  frequency  domain.  While  the  Fourier  Transform  S(u)  is  a  linear  function 
of  the  analyzed  signal,  the  Fourier  power  spectrum  |  S{u>)  |2  is  quadratic  to  the  signal  s(t). 
The  Fourier  Transform  S(u)  in  general  is  complex,  whereas  the  power  spectrum  |  S{u)  |2 
is  always  real.  The  Fourier  Transform  and  the  Fourier  power  spectrum  are  the  two  most 
important  tools  for  frequency  analysis. 

Before  investigating  the  case  of  an  accelerating  target,  it  is  helpful  to  review  the  case  of  a 
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non-accelerating  target. 

When  the  acceleration  is  zero,  and  the  noise  is  absent,  the  Fourier  power  spectrum 


S(u) 


_1_ 

N 


N- 1 


n~0 


2 


N 


N—l 

n=0 


2 


(9) 

(10) 


has  a  unique  maximum  which  is  attained  when  u  =  bx  mod  27t.  Assuming  — 7 r  <  <  71,  the 

maximizing  w  is  an  estimate  of  the  normalized  radial  velocity. 

Therefore,  the  target  can  be  detected  and  the  velocity  parameter  bx  can  be  estimated  by 
using  peak  detection  on  the  noisy  version 


r(w)  = 


N 


N-l 

J2  f{n)e~^n 


n—0 


2 


(ii) 


of  S(u),  that  is,  maximizing  r(oi)  with  respect  to  lo,  and  declaring  the  presence  of  a  target  if 
the  maximum  value  exceeds  a  threshold;  the  maximizing  u>,  which  is  now  a  random  variable 
due  to  the  presence  of  noise,  is  taken  as  the  estimate  of  bx.  This  method  is  equivalent  to  the 
GLRT. 


The  notion  of  output  SNR  of  the  method  comes  from  the  observation  that  under  the  ‘noise- 
only’  hypothesis  is  a  central  Chi-Squared  random  variable  of  two  degrees  of  free¬ 
dom,  for  all  u,  whereas  under  the  signal-and-noise’  hypothesis  )  is  a  non-central  Chi- 

Squared  random  variable  of  two  degrees  of  freedom  and  non-centrality  parameter  2 N  ^ 
[Whalen,  1971;  Meyer  and  Mayer,  1973]. 

We,  therefore,  consider 


(12) 

as  this  can  be  used  to  approximately  compute  the  probability  of  detection,  assuming  that 
the  maximum  of  T(o;)  occurs  at  u  =  bx.  Note  that  SNR0U^  is  a  monotonically  increasing 
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function  of  N.  More  specifically,  it  is  a  linearly  increasing  function  of  N.  Normalizing  this 
with  respect  to  SNRp^gg  gives 


SNR, 


gam 


^2BL,  (13) 

SNRpuise 

N.  (14) 
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4  The  performance  of  the  Fourier  method  in  detecting 
an  accelerating  target 

When  the  acceleration  is  non-zero,  and  the  noise  is  absent,  the  Fourier  power  spectrum  is 


SM  = 


1 

N-l 

2 

Jl 

N 

E  s(n)e-*~ 

71=0 

) 

17V— 1 

|2 

CL 

N 

eA(6i-^)n+5fc2n2) 

n=0 

Some  insight  into  this  function  can  be  gained  by  rearranging  it  as 


S(u)  = 


N 


n~  0 


where 


(15) 

(16) 


(17) 


bm  —  b  i  + 


N-l 


(18) 


is  the  normalized  median  velocity  over  the  time  interval  0  to  [N —  1  'jT  seconds.  For  simplicity 
of  discussion  assume  that  —  n  <  bm  <  n.  The  rearrangement  shows  that  S(u>)  is  symmetric 
with  respect  to  u  =  bm,  that  is,  S(bm  +  v)  =  S(bm  —  v),  and  the  spread  of  the  function  about 
u>  —  bm  depends  only  on  b2  and  N.  More  specifically,  when  considered  as  a  function  of  62, 
S(cu)  depends  only  on  the  magnitude  of  62. 

A  closed-form  expression  for  the  above  sum  seems  difficult  to  obtain.  But  numerical  eval¬ 
uations  suggest  that,  for  a  given  acceleration  b2,  there  is  a  threshold  value  Nt  of  N ,  below 
which  S(u>)  has  a  unique  maximum,  and  above  which  S(co)  has  at  least  two  equal  maxima. 
The  threshold  value  Nt  depends  on  the  acceleration  b2  alone.  Moreover,  Nt  is  a  decreasing 
function  of  b2.  Thus  Nt  implicitly  defines  the  maximum  allowable  acceleration  for  a  given 
N. 
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4.1  The  Region  N  <  Nt 


For  N  <  Nt,  by  symmetry  considerations,  the  unique  maximum  of  S(ui)  is  attained  when 
u  =  bm.  Thus 


max 

u> 


SM 


S(bm), 


a 

N 


Y^M»~^Y 

n— 0 


(19) 

(20) 


Therefore,  in  this  region,  we  have 


and 


SNR, 


SNRout  = 


SNRout 


SNRp^’ 
S(bm) 


cr 


1 

N 


Y^M-^Y 

n~  0 


(21) 

(22) 

(23) 

(24) 


Provided  that  SNRq^  is  adequate,  peak  detection  still  provides  a  means  of  detecting  an 
accelerating  target  and  estimating  its  median  velocity.  However,  as  expected,  SNR0U^  is  not 
a  linearly  increasing  function  of  N.  It  is  not  even  a  monotonically  increasing  function  of  N. 
Rather  it  is  a  concave  function  of  N  and  attains  a  maximum  value  in  the  region  N  <  Nt. 
The  optimum  value  of  N,  where  the  SNR0U^  is  maximum,  depends  on  the  parameter 
b2  alone.  Moreover,  Nopt  is  a  decreasing  function  of  62- 

For  N  <  N^t,  the  concentration  of  the  Fourier  power  spectrum  S(u)  about  u  =  bm  increases 
as  N  increases  towards  Nopt.  For  N  >  Nopt,  the  concentration  of  the  Fourier  power  spectrum 
S(u>)  about  u  —  bm  decreases  as  N  increases  towards  Nt. 


4.2  The  Region  N  >  Nt 

For  N  >  Nt,  the  Fourier  power  spectrum  S(u)  has  at  least  two  maxima  located  symmetrically 
about  u>  —  bm.  Moreover,  it  is  spread  symmetrically  about  u  =  bm,  and  the  amount  of  spread 
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increases  rapidly  with  N,  soon  covering  the  entire  interval  of  length  2n.  The  combination 
of  the  significantly  low  maxima  and  the  spectral  spreading  makes  it  practically  impossible 
to  discern  the  presence  of  the  target  or  estimate  its  spectral  location.  Therefore,  the  use  of 
peak  detection  as  a  method  of  detecting  an  accelerating  target  and  estimating  its  median 
velocity  is  precluded  in  this  region. 
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5  Numerical  Examples 


Figure  1  shows  how  the  Fourier  spectrum  changes  when  N  (referred  to  as  the  length)  is 
increased  while  the  acceleration  is  held  constant.  Here  we  can  see  the  spectrum  concentration 
first  increasing,  then  slowly  decreasing,  and  eventually  rapidly  decreasing  with  multiple  peaks 
appearing. 

Figure  2  shows  how  the  Fourier  spectrum  changes  when  the  acceleration  is  increased  while 
N  is  held  constant.  Here  we  can  see  the  spectrum  concentration  first  slowly  decreasing  and 
then  rapidly  decreasing  with  multiple  peaks  appearing. 

In  all  cases  the  median  velocity  was  held  at  zero  for  the  purpose  of  compact  representation, 
but  the  conclusions  are  independent  of  this  assumption. 
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6  Performance  Graphs 


Graphs  quantifying  the  performance  of  the  Fourier  method  in  detecting  an  accelerating  target 
are  given  in  this  section.  The  graphs  were  obtained  by  empirically  studying  the  Fourier  power 

spectrum  for  values  of  N  ranging  from  64  to  512  and  acceleration  62  ranging  from  0.00005 
to  0.0004. 

Figures  3  and  4  show  the  dependence  of  the  SNRgain  on  Length  N  of  the  signal  for  various 
values  of  the  acceleration  b2.  We  can  see  that  for  b2  >  0.00012,  the  failure  point  has  been 

achieved,  that  is,  Nt  <  512,  and  for  b2  >  0.00006,  the  optimum  point  has  been  achieved, 
that  is  Ngpt  <  512. 

Figure  5  shows  the  dependence  of  the  Optimum  SNRgain  on  acceleration.  Figure  6  shows 
the  dependence  of  the  Optimum  Length  Nopt  on  acceleration. 

Figure  7  shows  the  dependence  of  the  log  of  Optimum  SNRgain  on  the  log  of  acceleration. 
From  the  linear  least  squares  fit  we  can  infer  that 


10 l°Sio  SNRgain, opt  ~  5  log10  b2  +  4.537. 


(25) 


or 


SNR  •  ~  in04537A-0-5 

DiNngam,opt  ~  iU  b2  .  (26) 

Figure  8  shows  the  dependence  of  the  log  of  optimum  length  on  the  log  of  acceleration.  From 
the  linear  least  squares  fit  we  can  infer  that 


l°§io  N0pt  ~  0 . 5  logjg  b2  +  0.5606, 


(27) 


or 


N, 


opt 


10 


.0.5606  l -0.5 


(28) 
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Combining  these, 


SNRgainiopt  *  10 -°  luNvt.  (29) 

Thus  when  operating  at  the  optimum  point  N  =  the  loss  with  respect  to  the  non¬ 
accelerating  case  is  1.12dB.  This  is  also  the  loss  with  respect  to  optimum  coherent  integration. 

Every  effort  was  made  to  ensure  the  accuracy  of  the  above  results.  However,  caution  must 
be  exercised  in  applying  them  outside  the  studied  range,  especially  for  larger  accelerations. 

Figure  9  shows  the  dependence  of  the  log  of  threshold  SNRgajn  on  the  log  of  acceleration. 
From  the  linear  least  squares  fit  we  can  infer  that 


101og10  SNRgaill)0pt  »  -4.677 log10  b2  +  3.746,  (30) 

or 


SNRgain,0pt  «  100374V-4677.  (31) 

Figure  10  shows  the  dependence  of  the  log  of  threshold  length  on  the  log  of  acceleration. 
From  the  linear  least  squares  fit  we  can  infer  that 

log10  Nt  «  -0.5078  log10  fe2  +  0.6859,  (32) 

or 


Nt 


10 


0.6859l— 0.5078 
°2 


(33) 


Combining  these, 


SNR, 


gain, opt 


10 


—0.3745  jy0.9210 


(34) 
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7  Application  to  GLRT 


In  the  case  of  an  accelerating  target,  the  GLRT  method  seeks  to  maximize  A (Cl,c2)  = 


lEnJo1  r(n)e-^Cin+iC2"2) 

way  of  doing  this  is  to  evai 


N 


with  respect  to  ca  and  c2  [Whalen,  1971;  Kay,  1993].  A  crude 
iuate  A(ci ,  c2)  on  a  rectangular  grid  in  the  ci-c2  plane  and  choose 
the  maximum.  The  results  of  the  previous  section  imply  that  the  grid  spacing  in  the  c2 
dimension  cannot  be  less  than  a  critical  value  if  the  method  is  to  work  for  all  values  of  b2. 


To  see  this,  first  note  that,  for  a  fixed  c2,  maximizing  A(ci,  c2)  with  respect  to  Ci  is  equivalent 
to  using  the  Fourier  method  to  detect  a  target  with  initial  velocity  bx  and  acceleration  (b2-c2). 
Suppose  the  grid  spacing  in  the  c2  dimension  is  2 A.  Then  for  every  b2  there  is  a  c2  grid  value 
such  that  \b2  —  c2|  <  A.  For  this  fixed  c2  value,  the  noise-free  version  of  A(ci,  c2)  will  have 
a  unique  maximum  with  respect  to  ci  provided  the  threshold  length  Nt  corresponding  to 
acceleration  \b2  —  c2\  is  greater  than  N.  To  ensure  uniqueness  of  the  maximum  for  every  b2. 

we  require  A  <  A^j^,  where  A^.r^  is  the  acceleration  for  which  the  given  N  is  the  threshold 
length. 


The  SNRgam  is  at  most  N,  and  depends  on  |&2  —  c2|.  We  can  use  the  graphs  to  choose  A 
so  that  the  loss  of  SNRgajn  is  less  than  a  given  value. 
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8  Conclusions 


The  Fourier  transform  method,  or  Doppler  processing  method,  has  been  generally  used  in 
HF  radar  to  detect  targets  that  are  moving  with  constant  radial  acceleration.  Examples 
of  accelerating  targets  are  manoeuvring  aircrafts  and  missiles.  In  this  report  we  show  that 
there  are  limitations  and  shortcomings  to  the  Fourier  transform  method  to  detect  accelerating 
targets  because  of  the  phenomenon  known  as  Doppler  smearing. 

We  considered  using  the  Fourier  method  for  detecting  constantly  accelerating  targets  and 
estimating  their  median  velocities,  and  empirically  showed  that  there  are  decreasing  functions 
Nopt (62)  and  Nt{b2)  of  the  acceleration  b2,  termed  the  optimum  length  and  the  threshold 
length,  respectively,  with  the  following  properties: 

1.  in  the  region  N  <  Nt  the  Fourier  power  spectrum  of  the  noise-free  signal  has  a  unique 
maximum  at  the  median  velocity,  and  outside  of  which  it  has  at  least  two  maxima;  thus  Nt 
implicitly  defines  the  maximum  allowable  acceleration  for  a  given  N. 

2.  in  the  region  N  <  Nt  the  SNRgajn  varies  as  a  concave  function  of  N  achieving  a  maximum 
at  N  =  Nopt',  moreover  N^t  <  Nt. 

Formulae  were  given  for  Nopt  and  corresponding  SNR  in  terms  of  the  acceleration.  Similarly, 
Formulae  were  given  for  Nt  and  corresponding  SNR  in  terms  of  the  acceleration.  These 
formulae  might  lead  to  the  design  of  detectors  that  are  of  manageable  complexity  yet  have 
near  optimum  performance. 

As  an  application,  we  considered  a  GLRT  based  detector  that  applies  a  search  over  a  grid. 
We  showed  that  the  grid  spacing  must  be  less  than  a  critical  value.  We  also  showed  how  to 
choose  the  grid  spacing  so  as  to  achieve  a  given  SNRgajn. 
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